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Abstract 

This article concerns the dressing method for solving of multidimensional nonlinear 
Partial Differential Equations. In particular, we join hierarchy of matrix Burgers type 
equation with hierarchies of equations integrable by the Inverse Spectral Transform (1ST). 
Example of resonance interaction of wave packets in (3+l)-dimensions is given. 



O 1 Introduction 

: 

o ■ 

It is well known that different versions of the dressing method [1, 2, 3, 4, 5] are very successful 
tools for solving of so called completely integrable nonlinear Partial Differential Equations 
(PDE). These equations, in turn, have wide application area in different branches of physics, 
such as hydrodynamics, plasma physics, superconductivity, nonlinear optics. Although till 
recently dressing methods have been used only for PDE integrable by the Inverse Spectral 
Transform (1ST), it has been shown in [6] that there is another type of equations (maybe not 
integrable by 1ST) which admit properly modified dressing procedure for construction of large 
manifold of their solutions. But the technique proposed there left opened many questions. For 
instance, it is not clear whether derived nonlinear PDE can be linearized by some substitution 
[7, 8]. Also it was difficult to characterize the manifold of available solutions. 

In this paper we replace algebraic operator with integral one, generalize system of equations 
introducing the set of additional parameters (independent variables of nonlinear PDE) and 
modify significantly the algorithm given in [6]. This allows us to simplify description of PDE's 
properties, exhibit more information about solution manifold as well as relations to the classical 
solvable (both linearizable and integrable by 1ST) PDE. Although all statements following 
hereafter can be proved, we omit most of the proofs for the sake of brevity. They will be given 
in different paper. 

Thus, the basic object is the following N x N matrix integral equation 

$ = $(A, fi;t) = J #(A, v\ t)U{v, //; t)dv = * U, (1) 
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where A = (Ai, . . . , \n x ), = (a*i? ■ ■ ■ ■> f^N^), v = (vi, . . . , vn v ) are vector variables with different 
length in general; integration is over whole space D u of the appropriate vector parameter; $, 
\I/ and C/ are N x N matrix functions of arguments. Star means integration over space of 
inner variable: f * g = f f(\\, . . . , v)g{y, . . . , fj, n )di>. We require that \& is invertible operator, 

i.e. equation (1) can be solved uniquely for the function U. By definition, operator A(X, (J,) 

is invertible, if there are operators A L 1 (X, fx) and Aj^iX, fi) such that J A(X, v)A~£{v,n)dv = 

n, 

J ^4£ 1 (A, v)A{y, n)dv = S(X—/j,). The functions <3> and \l/ are related by means of the compatible 

D v 

system of linear integral-differential equations, which, on the other hand, introduce (infinite) 
set of additional parameters t — (ti,t 2 ■ ■ ■ ) (independent variables of nonlinear PDE): 

Mi *^ = J2 L ik*®*Cki, C kl = C ki (X,fi-t), i = l,2,... (2) 

k 

where Mj = M,(A, /i; d tl , d t2 , . . . ) are first order and Lj k = Lj k (X, /i; d tl , d t2 , . . . ) are arbitrary 
order linear differential operators with matrix coefficients depending on A and This overde- 
termined system together with its compatibility condition defines \1/ and $. Finally, the same 
compatibility condition with substitution $ from the eq.(l) results in nonlinear PDE whose 
solution is expressed in terms of U . 

After this preliminary discussion we derive some general equations using the following sim- 
plified version of the system (2): 

m u = S % *^*C l , i = 1,2,..., (3) 

where Si(X,fi;t) and Ci(X,fi;t) are known functions of t, which will be seen later. The com- 
patibility condition for the system (3) has the form 

S it . * $ * Ci - S h . *<f>*Cj + Si*$* C it . - Sj * $ * C ju + (4) 

St * $ tj . *Ci- Sj * $ t . * Cj = 0, 

which is linear system of compatible integral-differential equations for the function $. Solving 
this equation, substituting result in (3) and integrating it, we obtain the expression for ^: 
^(A, pi; t) = d^ 1 1 (S 1 *^*C 1 )(X, /i; t) + E(X, /jl) + F(X, /i; t 2 , t 3 , . . . ). Here E is invertible operator, 
function F provides compatibility of the system (3). Being invertible, operator \l/ provides 
unique solution to the eq.(l). 

On the other hand, eq.(4) may be given another form after substitution eq.(l) for $ and 
(3) for ^ 

S it . * * * U * d - S ju * * * U * Cj + Si * * * U * C it . - Sj * * * U * C jt + (5) 
S i *(S j *V*U*C j *U + ty* U tj ) *Ci- Sj*(Si*V *U *d*U + * *U U )* Cj = 0, 

which is nonlocal equation quadratic in U . It may result in nonlinear PDE for the dependent 
variables, expressed in terms of U , St and C{. To provide this possibility we must impose 
specific dependence of the functions Si and Cj on their arguments, for instance, in accordance 
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to the following set of relations: 

S t (\,fi;t) = S(\,fi), $ = 1 S*$ + x, X = x(\Mt), (6) 

Ci(\,n;t) = J A i (X,u)p 1 (u;t)p 2 (fj)du + c 1 (X)B i c 2 (fJ>;t) = A t *pi{t)p 2 + c 1 B i c 2 (t), (7) 
d v 

Ai*d = ciB h (8) 

where A is invertible operator, Ai * Aj = Aj * A { , [B h Bj] = and [*, *] means commutator of 
two matrices. Eqs.(6-8) split eq.(4) into the following set of three integral-differential equations 
for $, pi and c 2 : 

S * $ t . * Ai - S * $ tl * Aj = 0, (9) 
A x * ( Pl ) tj - Aj * (px) tl = 0, (10) 
fi lC2t3 - B jC2tl =0. (11) 

Then one has the following nonlinear equation instead of (5) 

ty*(U tj *Ai-U tl *A j + U*C j *U*A 1 -U*C 1 *U* Aj)+ (12) 
X tl * Aj -Xt ] *A 1 + x *(C l *U* Aj -Cj*U* A x ) = 0. 

Note, that reduction leading to the equations introduced in [6] will be discussed in different 
paper. Here we consider two other examples of multidimensional systems. First of them 
(Sec. 2.) represents combination of linearizable (Burgers type) and completely integrable (n- 
wave) (3 + l)-dimensional systems, having solutions depending on arbitrary functions of three 
variables. Second example (Sec. 3) is another generalization of the matrix n-wave system [9]. 
Properly introduced multiple scales expansion of this system results in the multidimensional 
((3+l)-dimensional in our case) equation describing resonance interaction of wave packets. Its 
solutions may depend on arbitrary functions of two variables. Both examples have extension 
into (n + l)-dimensions with arbitrary n. 



2 Generalized hierarchy of linearizable (Burgers type) 
and integrable by 1ST (n-wave) systems 

In this section S(\,p) = 5(X — //), x — 0> Aj = A* ■ • ■* A = A>, Bj = B j , where A(X, yi) is 

3 

invertible operator and B is nondegenerate constant matrix. Thus ^ = $ * d with given 
by (7). After applying operator from the left to the eq.(12) one results in 

Ej =U tj * A + U * A j * pip 2 *U * A + U * ciB j c 2 *U * A- (13) 
(U tl * A j + U * A * pip 2 *U * A j + U * ciBc 2 *U * A j ) = 

We may derive nonlinear system for the functions 

u = p 2 *U * ci, q n = p 2 * U * A n * p ± , v n = d n c 2 *U * c ± , w nm = d n c 2 *U * A m * p ± , (14) 
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which has the following "short" form 

p 2 * Ej * ci = 0, p 2 * Ej * A n * Pl = 0, <9 n c 2 * Ej * ci = 0, <9 n c 2 * * A m * Pl = 0, (15) 
or, extended form 

u t .-u tl B j - 1 + q j u-q 1 uB j ~ 1 + uB j v - uBvqB^ 1 = 0, (16) 
q nt . - q n +j-i tl + qjqn - qiQn+j-i + uB 3 w 0n - M5w 0(n+i _i) = 0, 
v ntj ~ Vn^B 3 ' 1 + w nj u - WniUB 3 ' 1 - B 3 ~ 1 v n+1 + v n+1 B 3 ~ l + (17) 

VnBivo-VnBvoB 3 - 1 = 0, 

W mntj - W m (n+j-l) tl + W mj q n - W m iq n+j -i - B 3 ~ 1 W(r n +l)n+ (18) 
W( m+ l)( n+J -_l) + v m B 3 w 0n - V m Bw ( n +j-l) = 0. (19) 

The complete system of pure PDE is represented by the following set: eq. (16) with j = 2, eqs. 
(17) and (18) with j = 2, 3, eq. (19) with j = 2, 3, 4. Thus, this system is (3+1) -dimensional. 
It may be given the compact form if one introduces column of matrices u,q n ,v n ,w nm : x — 
[u, q 1 ,q 2 ,...,v 1 ,v 2 ,..., w 00 , w 10 , w 01 ,... ] T : 

4 

^ ^ ^ ^ VlijmndtiXmn ~t~ ^ ^ T^ijklmnX.klX.mn 0, (20) 
1=1 mn klmn 

where Vuj mn and T^i mn are constants expressed in terms of the elements of the matrix B. 

Physical application of the eqs. (16-19) is not found yet. In particular, it reduces into the 
following (2+l)-dimensional systems: 

1. Matrix Burgers type system (i.e. linearizable) for the function q , if c\ = or N = 1. 

2. Matrix n-wave equation (n = N(N — l)/2) for the function v , if pi = 0. 

2.1 Construction of solutions 

First, one needs to solve the system (8-10) for the functions ci, $, p±, c 2 : 

A*ci = dB, (21) 
$(A,a«;0 = J J ^o(\^k)e m{u ' At) (f)o(^^k)dkdu, (22) 

Pi(X;t) = J j 'p (X,iy-k)e^ k ^p 10 (iy-k)dkdiy, (23) 

c 2 (A;0 = J e k ^ BHi c 20 (X;k)dk (24) 

where rji(fj,] k,t) = X^=i k)tj, i — 1,2, [%-,?7ifc] = 0, det^) 7^ 0. Parameter k is complex 
in general, integration is over whole complex plane of this parameter. Function c 2 q is 
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arbitrary and functions 0o and p solve the following system: 

77ij(i/;A;)0o(^A*;^) = j Vi(j-i)(v, fyfaiy, v i\ k)A(i/ u fj^dvx, (25) 

p (\,is;k)r} 2j (is;k) = J A(X, v^pofa, v\ k)r) 2 {j-\){v\ k)dv x , j = 2, 3, . . . . (26) 

Functions 77^1 (j = 1,2) are arbitrary, while r]j n with n > 1 provide compatibility of eqs. (25) 
and (26). 

Now one can integrate (3) to get \& (j = 1): 

tf(t) = E + di 1 1 [$(t)*A*p 1 (t)p 2 + $(t)*c 1 c 2 (t)}, (27) 

where .E = E(X,/i) is invertible operator independent on t. For instance, E{X,ji) = 5(X — fi). 
Next, find U from (1): U = * $. In general, operator can be constructed only 
numerically, unless <E> is degenerate ($ (A, /i; A;) = $ 01 (A)<l ) o2(/ i ; ^)) and explicit form for 
Ej^ 1 is known. In this case ^l/ -1 may be found analytically, following the procedure proposed, for 
instance, in [10], where <9-problem with degenerate kernel has been solved. Similarly, eqs. (25) 
and (26) can be solved numerically, unless A has the following structure: A(X, /i) = A (X, //) + 
Aji(X)Aj 2 (fj,), where operator A is invertible with known analytical form for Aq]^ . For 
instance, A (X, //) = 5(X — //). Then compatibility condition of the system (25) and (26) 
produces dispersion relations in the form r)i n (n; k) = 7711 (/i; k) i*i[(0o * k), % — 1, 2, . . .], 

^2n(A*; fc) = ^2[(^-i2 * Po)(a*; fc), * — 1? 2, ... ] r) 2 i(ii;k), where Fj are given matrix functions of 
matrix arguments. 

Finally, one can show that solutions of our (3+l)-dimensional system (16-19) constructed 
in accordance with definitions (14) may depend on arbitrary functions of three real parameters, 
for instance, t 1: t 2 and t 3 . This is owing to the factor ^ * A * p x . 

3 Resonance wave interaction in (3+l)-dimensions 

In this section we consider eqs. (6-11) with S(X,/i) ^ 5(X — /i), x 7^ 0, p\ = and \J/ ti = S , *$*C i . 
It is convenient to apply operator c\ to the eqs.(l) and (9) from the right, giving them the form: 

S*§ + x = \&*U, U = U*Ci, x = x*ci, (28) 
(S*Q) t . = (S*Q) tl Bj, B 1 = I, (29) 

/ is identity matrix, Bi are diagonal matrices and iV > 4. Let Xtj(X;t) = x(X;t)dj, where dj 
are constant matrices. We will need the following notations: bj = a±Bj — a,j , V$ = c 2 * U and 
Vi = c 2t! * U. Nonlinear eq.(12) gets the following form after applying operator c\ from the 
right: 

* (d tj U - d tl UB 3 + U[Bj, Vo\) + x(bj ~ [Bj, V ]) = 0. (30) 

Now assume that det(6j — [Bj, Vq[) 7^ for all j and use two equations (30) with indexes j and 
k, j 7^ k to eliminate function x- After applying operator c 2 * from the left to the resulting 
equation, we receive: 

(d tk V - d tl V B k + [V ± , B k ] + V [B k , V ])(b k - [B k , Vo})' 1 = (31) 
(d tj V - d tl V Bj + [V u Bj] + V [Bj, V })(bj - [Bj, Vo})' 1 . 
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Next, let us introduce different scales for variables t k) V , V\. d tk — > ed tk , V — ev, V\ — e 2 v\. 
Keeping only leading terms, we get from the eq.(31): 

E k = v tl (Bjbj 1 - B k b- V ) + v tk b- 1 - vtjbj 1 + (32) 
[vi, B k ]b? - [ Vl , B^bj 1 + v[v, B^bj 1 - v[v, B k ]b? = 0. 

Thus the complete system is represented by the pair of equations (32), E k and E n , k ^ n. One 
can see that the following combination of these equations has no function v\ and contains only 
off-diagonal elements of v. 

E k (B n b~ l - Bjbj 1 ) - E n (B k b k l - B.bj 1 ) + (33) 

Bj(E k — E n )b~j 1 — B n E k b~^ + BkEJ)^ 1 = 0. 

Let j = 2, k = 3, n = 4 and write this equation in the following form 

4 

^2 S nijdt n Vij + ^ T ikjVikV k j = 0, i ^ j, (34) 
71=1 k:kj^i^j 

where s ki j, and T ik j are constants, expressed in terms of the elements of the matrices Bj and 
by If Vij are real, then this equation describes resonance interaction of wave packets. 

Reduction t k — > it k , = v jh with real s nij and T ikj , s nij = s nji , T ikj = T jki (bar means 
complex conjugated value) transforms the (3+l)-dimensional eq.(34) into (2+l)-dimensional 
n-wave equation with independent variables r k — t k + h, k — 2, 3, 4. 

3.1 Construction of solutions 

In this section we give the algorithm for construction the solution Vo to the eq.(31). 
Solutions of the eq.(29) and expression for x have the form: 

S*$(A) = J $ (\,k)e k ^ Bntn dk, c 2 (A) = J e k ^" Bntn c 20 (\,k)dk, (35) 

x(A)= X o(A)e^ a " t ". (36) 
To find \1/ we integrate eq.(3) (j = 1, remember that Bi = I): 

*(A,a*) = | J M\,k)e^^ B ^c 20 (^q)^ + 5(\- fi). (37) 

Thus 

U(\)=S*$(\)- J y$ o (A,A;)e^ B ^ (fc+9) 0(g)^ + x(A), = c 20 * £/ (38) 

and 

Vb = ca * U. (39) 

Unknown function 0, related with U, can be found only numerically in general case, unless 
functions $(A, k) is degenerate [10]. Eq.(39) shows that V may depend on N x N matrix 
function of two real variables, for instance, t\ and t 2 . 

Regarding the multi-scale expansion given by eqs.(32), one should replace t n — > et„, in 
formulae (35-39) and take arbitrary functions <3>o and xo proportional to e. Thus Vq ~ e. 
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4 Conclusions 



Working with dressing methods we underline two directions: (a) increase of dimension of solv- 
able nonlinear PDE and (b) provide rich class of their solutions. Nonlinear PDE derived with 
our algorithm admit infinite set of commuting flows corresponding to different parameters tj. 
Since general equations are ruther complicated (see (16-19), (31)), the reasonable problem is 
construction of their reductions, which would exhibit physical application of these systems. 
Another way is multi-scale expansion of general systems, which in our case reveals (3+1)- 
dimensional equation describing resonance interaction of wave packets (see eq.(32) and (34)). 

Author thanks Prof. S.V.Manakov for discussion. This work was supported by RFBR 
grants 03-01-06122 and 1716.2003.1. 
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